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I,  INTRODUCTION 


When  the  geometrical  theory  of  diffraction  (GTD)  [1,2]  is 
employed  to  analyze  the  problem  of  high  frequency  plane  wave 
scattering  by  smooth  convex  surfaces,  a  modification  of  this 
method  is  required  for  analyzing  the  fields  in  the  vicinity  of 
the  shadow  boundaries  associated  with  this  problem.  This  modi¬ 
fication  of  the  GTD  forms  the  subject  of  the  present  paper. 


In  the  GTD  analysis  of  the  diffraction  of  waves  by  a  smooth 
convex  cylinder  [3],  the  total  field  exterior  to  it  is  associated 
in  a  simple  manner  with  the  usual  incident  and  reflected  rays  of 
geometrical  optics  together  with  the  surface  diffracted  rays 
introduced  by  Keller.  The  geometric  optical  rays  are  present 
only  in  the  lit  (or  illuminated)  region,  so  that  the  surface  dif¬ 
fracted  rays  alone  account  for  the  field  in  the  shadow  region 
(these  surface  diffracted  rays  are  also  present  in  the  illumin¬ 
ated  region  if  the  cylinder  cross-section  is  closed).  However,  pure 
ray  optical  field  descriptions  fail  at  and  near  shadow  boundaries, 
and  at  caustics  of  ray  systems;  consequently,  the  GTD  ray  solution 
for  the  problem  illustrated  in  Fig.  1  is  valid  in  regions  I  and  III, 
but  it  fails  in  regions  II,  IV,  V,  and  VI.  Region  II  constitutes  the 
penumbral  or  transition  region  adjacent  to  the  shadow  boundary  (SB). 
The  angular  extent  of  this  transition  region  is 


0 
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where  k  refers  to  the  wavenumber  of  the 


surrounding  medium  which  is  assumed  here  to  be  free  space,  and 
Pq(Ql)  is  the  radius  of  curvature  of  the  cylinder  at  the  point  of 
grazing  incidence,  i.e.,  at  Qi.  Regions  IV,  V,  and  VI  are  usually 
referred  to  as  the  surface  or  caustic  boundary  layer  regions.  In 
shadow  regions  IV  and  V,  the  surface  of  the  obstacle  is  a  caustic 
of  the  surface  diffracted  rays;  whereas,  region  VI  is  in  the  neigh¬ 
borhood  of  Qi  which  is  a  caustic  of  the  reflected  ray  for  grazing 
incidence.  The  mathematical  representations  for  the  field  within 
the  transition  regions,  the  surface  boundary  layer  regions,  and 
also  their  common  regions  of  overlap  are  complicated  by  the  fact 
that  they  must  change  rapidly  but  smoothly  from  one  form  to  an¬ 
other  across  these  regions. 


It  is  evident  that  the  usefulness  and  versatility  of  the  GTD 
solution  is  further  enhanced  if  it  can  be  modified  within  the  region 
where  it  fails;  such  modifications  may  be  achieved  by  introducing 
uniform  solutions  which  not  only  remain  valid  in  one  or  more  of 
these  regions,  but  which  reduce  to  the  conventional  GTD  results 
exterior  to  these  regions. 

The  development  of  a  simple,  uniform  solution  which  remains 
valid  in  region  II  that  is  exter'or  to  the  urface  boundary  layer 
is  considered  in  this  paper.  Th  ‘-’nrk  '  iiotivated  by  the  need 
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Fig.  1.  Rays  associated  with  the  plane  wave  scattering 
by  a  smooth  convex  cylinder. 
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to  analyze  the  radiation  patterns  of  antennas  mounted  on  the  wings 
or  tails  of  aircraft  [4,5].  The  radiation  patterns  of  such  antennas 
are  significantly  affected  due  to  shadowing  by  the  aircraft  fuselage 
surface  which  is  generally  smooth  and  convex  shaped.  In  many  such 
wing  or  tail  mounted  aircraft  antenna  configurations,  a  knowledge 
of  the  field  in  the  transition  (or  penumbra!)  region  is  essential 
for  the  purposes  of  antenna  design.  One  can  develop  a  solution  for 
the  transition  region  via  an  asymptotic  high  frequency  analysis  of 
a  canonical  problem  which  closely  models  the  local  geometry  of  the 
convex  cylinder  at  . 

In  this  paper,  the  canonical  problem  of  plane  wave  scattering 
by  a  perfectly-conducting  circular  cylinder  is  solved  by  introducing 
certain  asymptotic  approximations,  which  are  based  on  physical 
considerations,  to  obtain  a  uniform  asymptotic  result  that  is  valid 
in  the  transition  region  exterior  to  the  surface  boundary  layer. 

This  uniform  result  yields  a  finite,  continuous  solution  for  the 
total  field  across  the  SB  within  the  transition  region.  Exterior 
to  this  transition  region,  it  reduces  to  the  usual  GTD  solution 
in  terms  of  surface  diffracted  and  geometric  optical  ray  fields 
for  the  shadow  and  lit  zones,  respectively.  Furthermore,  this 
result  is  in  a  form  which  is  simple  and  convenient  for  engineering 
applications.  Upon  invoking  the  principle  of  the  local  character 
of  the  high  frequency  fields  in  the  penumbra,  one  can  systematic¬ 
ally  modify  the  uniform  result  for  the  circular  cylinder  to  treat 
the  electromagnetic  scattering  by  arbitrary,  smooth  convex  surfaces 
for  both  two-  and  three-dimensional  problems.  These  extensions 
will  be  presented  in  detail  in  a  separate  paper. 

The  analysis  of  the  canonical  problem  is  indicated  in  section 
II.  Numerical  results  based  on  this  uniform-GTD  analysis  are  com¬ 
pared  in  section  III  against  the  results  obtained  from  an  exact 
eigenfunction  series  solution  for  the  same  canonical  problem  where, 
it  is  seen  that  the  present  uniform  result  is  accurate  not  only 
within,  but  also  exterior  to  the  transition  region.  A  surmary  and 
discussion  follows  in  section IV.  Before  proceeding  with  the  analysis, 
the  present  work  will  be  compared  with  the  work  of  others  on  this  subject. 

The  problem  of  asymptotically  estimating  the  fields  within  the 
transition  regions  associated  with  the  diffraction  of  waves  by 
smooth,  convex  surfaces  has  had  long  standing  interest,  especially 
in  connection  with  the  theory  of  radio  wave  propagation  around  the 
earth.  The  fundamental  work  of  Fock  on  this  subject  [6]  appears  to 
be  the  most  significant;  in  particular,  his  analysis  of  the  fields 
near  a  shadow  boundary  and  within  the  Fresnel  region  of  a  diffracting 
spherical  surface  [7]  is  of  direct  relevance  to  the  problem  con¬ 
sidered  here.  Other  notable  works  which  are  also  directly  related 
to  the  present  analysis  are  due  to  Wait  and  Conda  [8],  Logan  and 
Yee  [9,10],  and  Ivanov  [11].  An  extensive  bibliography  on  various 
aspects  of  the  problem  of  diffraction  by  a  smooth,  convex  surface 
may  be  found  in  the  papers  by  Logan  [9,10],  and  Borovikov  and 
Kinber  [12].  Indeed,  the  present  work  draws  upon  the  earlier, 
important  contributions  of  Fock  and  Logan. 
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In  Fock's  analysis  [7],  the  wave  equation  within  the  transition 
region  is  solved  approximately  via  the  method  of  parabolic  equations; 
this  solution  is  initially  expressed  in  terms  of  a  canonical  integral 
introduced  earlier  by  Fock  [13a, b].  For  observation  points  near 
the  shadow  boundary  and  the  diffracting  surface  ^but  outside  the 
surface  boundary  layer),  this  canonical  integral  is  reduced  by  Fock 
to  integrals  which  are  functions  of  only  a  single  parameter;  in 
particular,  the  solution  is  expressed  in  terms  of  a  Fresnel  integral 
and  a  Pekeris-type  integral  [9,10],  both  of  which  are  well  tabulated. 

Wait  and  Conda  [8]  analyze  the  field  in  the  transition  region 
scattered  by  a  circular  cylinder  via  an  asymptotic  evaluation  of 
the  radiation  integral  which  employs  the  well-known  Fock  currents 
[13a]  for  this  problem;  their  asymptotic  solution  near  the  shadow 
boundary  is  also  given  in  terms  of  Fresnel  and  Pekeris  integrals 
as  done  earlier  by  Fock  [7].  In  addition,  the  result  given  in  [8] 
is  valid  not  only  near  the  surface  (outside  the  surface  boundary 
layer)  as  in  [7],  but  also  in  the  far  zone.  The  analysis  in  [8]  is 
an  extension  and  generalization  of  the  work  of  Goriainov  [14]. 

The  results  in  [7],  [8],  and  [14]  are  quite  accurate  in  the 
illuminated  region  very  close  to  the  SB,  but  in  general  they  do  not 
join  smoothly  with  the  geometrical  optics  field;  whereas,  the  result 
presented  here  does  reduce  to  the  geometrical  optics  field  in  the 
illuminated  region.  Also,  the  present  result  is  valid  in  both  the 
near  and  far  zone  as  in  [8];  furthermore,  exactly  on  the  shadow 
boundary,  it  agrees  with  the  results  of  [7],  [8],  and  [14]  which 
are  accurate  near  the  shadow  boundary.  Logan  and  Yee  [10]  give  a 
result  for  the  illuminated  part  of  the  transition  region  in  terms  of 
the  canonical  integral  of  Fock  such  that  it  reduces  to  ihe  geometrical 
optics  field  far  from  the  shadow  boundary.  This  result  appears  to 
have  been  constructed  by  redefining  certain  parameters  in  Fock's 
canonical  integral  on  an  ad  hoc  basis. 

Ivanov  [11]  analyzes  the  same  canonical  problem  as  in  this  paper 
via  the  method  of  parabolic  equations;  his  solution  which  is  valid 
in  the  transition  and  the  surface  boundary  layer  regions,  and  which 
reduces  to  the  GTD  solution  outside  the  transition  region,  is  also 
given  in  terms  of  an  integral  related  to  the  canonical  integral  of 
Fock.  The  canonical  Fock  integral-type  solutions  [10,11]  are 
complicated  for  numerical  evaluation  in  the  transition  region  except 
when  the  field  point  is  either  very  close  to  the  shadow  boundary, 
or  is  on  the  surface;  in  these  exceptional  cases,  they  reduce  to 
well  tabulated  integrals  which  are  functions  of  only  a  single 
parameter. 

On  the  other  hand,  using  a  somewhat  less  formal  analytical 
procedure  than  the  parabolic  equation  method  employed  in  [7], 

[10],  and  [11],  the  present  solution  is  given  in  terms  of  the 
well  tabulated  single  parameter  Fresnel  and  Pekeris  integrals  for 
the  entire  transition  region  (exterior  to  the  surface  boundary 
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layer).  In  this  sense,  the  present  result  Indicates  a  simplifica¬ 
tion  over  the  results  in  [7],  [10],  and  [11],  respectively,  for 
field  points  outside  the  surface  boundary  layer;  furthermore,  this 
simplification  is  achieved  without  any  apparent  loss  of  accuracy  as 
confirmed  by  the  numerical  calculations  in  section  III. 


II.  ANALYTICAL  DETAILS 

An  asymptotic  high  frequency  analysis  of  the  canonical 
problem  of  electromagnetic  plane  wave  scattering  by  a  perfectly- 
conducting  circular  cylinder  in  free  space  is  presented  in  this 
section;  both,  the  TE  and  the  TM  polarization  cases  are  considered. 
The  geometrical  configuration  of  this  canonical  problem  is  illus¬ 
trated  in  Fig.  2. 


Fig.  2.  Geometry  of  the  canonical  problem. 


Let  the  incident  plane  wave  field  be  denoted  by  y\  and  the 
scattered  field  by  1)5;  then,  the  total  field  is  U  =  Uj  +  exterior 
to  the  cylinder.  The  total  field  U  corresponds  to  a  z-directed 
electric  field  if  the  excitation  is  TM^  polarized;  whereas,  it 
corresponds  to  a  z-directed  magnetic  field  if  the  excitation  is  TE^ 
polarized.  The  form  of  the  incident  field  U''  is: 


n) 


eikx 


The  total  field,  U  satisfies  the  reduced  wave  equation 
(2)  +  k^)  U  =  0  , 
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where  v  is  the  usual  2-D  Laplacian  operator  in  cvlindriral 
coordinates  (p,(p).  The  scattered  field  l)S  satisfies  iho  Somnioi- 
feld  radiation  condition  for  an  time  dependence  (which  is 

assumed  and  suppressed),  and  the  total  field  U(p,(;))  satisfies  the 
following  boundary  conditions  on  the  cylinder: 


(3) 


Q 


U(p  ,(|) ) 

=  0  ,  at  p  =  a;Q  =  J 

-  -J 

1 

1 ,  for  the  TM^  case 

3 

for  the  TE^  case. 


The  solution  to  (2)  subject  to  the  radiation  condition  and  (3) 
is  [15] 


(4a) 


or 


U  = 


a>-ie  r 

■i  j  dv 
-“-ie 


Q  J„(k3)  foN 


,  cos  V  (tt  -  U  I  ) 
Sin  VJT 


IV 


(4b) 


OO-I  C 


U  =  -^ 


dv 


-«-le 


H  (1) 


(kp)  - 


Q  H^^^Va) 


fp) 


COS  V  (ir  -Ul) 


Sin  viT 


with  e^O,  and  e-*^0.  The  quantities  Jy,  ^  are  the  usual 

cylinder  Bessel  and  Hankel  functions.  It  is  noted  that  lm(v)<0  on 
the  contour  of  integration  in  (4a)  and  (4b). 

The  above  integral  expressions  for  U  may  be  rewritten  by 
employing  the  decomposition 


(5)  cos  v(7T  -  Ul)  =  cos  V(p  +  i  e"’''l'^l  sin  vt 


Thus,  for  large  ka  (i.e.,  ka>>l),  one  obtains 


(6a) 


U  % 


“-Ie 


•lE 


dv 


Jv(kp)  - 


Q  dv(ka) 

Q  H/^^(ka) 


(2) 


(kp) 


■iv[ki4] 


6 


ro|=4 


or 


(6b) 


“  '\i  2 


I 


dv 


Q  H^^^^(ka) 


-ivEUI 

e 


with  0<|4ij<7r  .  In  going  from  (4)  to  (6)  via  (5),  the  terms  containing 
e''''^  cos  v<(>  are  discarded  because  they  can  be  shown  to  correspond 
to  the  multiply  encircling  GTD  surface  ray  modes  which  contribute 
negligibly  for  large  ka.  Within  the  transition  region,  v  0  (ka) 
so  that  J  (ka)  and  H  become  for  large  ka: 


(7a) 


(7b) 


Jy(l<a)  %-J- V{m"^[v-ka])  ; 


(1) 

(2), 


H^UJ(ka)  W^(m'’[v-ka]) 

m/rr  2 

il/3 


fka  \ 

where  m=(^  I  ,  and  the  quantities  V  and  are  the  Fock-type 

^  ^  2 
Airy  functions*  [9].  Introducing  the  change  of  variable 
T  =  m"^[v-ka]  now  allows  one  to  express  (6)  in  terms  of  (7)  as: 


(8a) 


U  m 


e 


-Ott-l  j; 


-iv(T)'l' 


Tt-r/3 

e  ;  2  iV(t)  = 


1.2 

.  2it 


runs  from  “>e  to  “-ie  ;  whereas,  r2 


W-|(t)  -  W2(t).  The  contour 

is  the  complex  conjugate  of 


7 


or 


“^1  e 


(8b) 


?  1 


dT 


_<«ui  e 


Q  W,  ( ^) 


■iv(T)i|j 


with 

(8c) 


Q 


for  the  TM^  case 
for  the  TE^  case 


Also,  v(t)  -  ka  +  m  t,  =  k  1-  -5  ,  and  Iin(T)<0  in  the  integrands  of 
(8a)  and  (8b). 


The  integral  representations  for  U  in  (8)  will  next  be 
asymptotically  approximated  for  the  shadowed  zone  of  the  transi¬ 
tion  region  in  part  A  of  this  section;  whereas,  the  asymptotic 
approximation  for  U  in  the  illuminated  zone  of  the  transition 
region  will  follow  in  part  B. 


A.  Transition  Region  -  Shadow  Zone 

One  begins  by  rewriting  (8)  in  a  manner  analogous  to  that 
done  by  Fock  [7]  in  the  decomposition  of  a  canonical  integral 
representing  the  solution  of  his  parabolic  equation  for  the 
problem  of  diffraction  by  a  sphere;  namely,  U  is  written  as  a 
superposition  of  two  integrals,  one  being  the  integral  in  (8a) 
over  the  contour  from  0-ie  to  --ie  (defined  as  r:.r.lour  Co),  and 
the  other  being  the  integral  in  (8b)  over  the  contour  from 
-»-ie  to  0-ie  (defined  as  contour  C-]).  Thus, 


(9)  U  +  I2 

where 


(10a) 


In  = 


C1+C2 


Jv(^)^kp) 


,-1v(t)i</ 


di  - 


m 

1 


,-iv(T).|y 


'1 


dt 


8 


for  the  transition  region  adjacent  to  SB^  i 
is  similar. 

The  first  integral  in  I-i  can  be  solved 
integral  representation  for  Jy(kp),  namely. 


o,(kp)  =  ^  j 


and  interchanging  the  orders  of  integration 


(11)  m  j  d,  = 


ikp  cos  4i  X  . 

,  for  0<<)><iT  and 


Clearly,  (11)  is  just  the  incident  field, 
also  made  of  the  following. 


m 

where  6(i/;-6)  is 


.  .  1  j  dT  e-'"'*-®’’ 

C  +C 
^1  ^2 

the  Dirac  delta  function. 


Next,  the  second  term  on  the  R.H.S.  of 
expressed  in  terms  of  the  standard  integral 
Hviljdip)  as 

-7  (  H  W(k»)  d, 

=  -S  J  d,  }  dB 

Cl  Cj 


I  the  domain 
by  employing  the  usual 

1  +  iv0 

to  yield 

^-ikp  sin  4< 
e 


n  deriving  (11 )  use  is 


» 


(10a)  may  be  re¬ 
representation  for 


g-ikp  sin  B  +  i v(t )e 


which  for  Im  (t)<0  reduces  to 


(12) 

with 


-?|  H,<2)(kp)  d,  . 

^1 

'i‘  =  ‘i‘  -  ^  t  I'/'hj  for  0<(|><ir. 


g-ika  ^ 
2  ir  i 


dg 


-ikp  [sin  s  - 

C  p 

-  3 


i 


Fig.  4.  Contours  of  integration 
for  Bessel  functions. 


The  contours  of  integration  and  Cg  occurring  in  the  defini¬ 


tions  of  Jy(kp)  and 


(2) 


(kp)  are  Indicated  in  Fig.  4,  The  integrand 
on  the  R.H.S  of  (12)  has  a  saddle  point  at  b  =  where  B^  =  cos"^ 


Since  B^  is  a  constant  for  fixed  values  of  p  and  a,  it  is  conveniently 


determined  by  setting  (p  =  <P^  (for  which  the  observation  point  is  on 


SB^);  consequently,  B^  =  7  ~  '^s  ~  '*’s  "  7  ’  shown  in  Fig.  2.  Also, 
the  integrand  in  (12)  has  a  pole  at  B  =  Bp,  where  Bp  =  ^  =  ii>  -  . 


Let  0  H  ~  then 


(13)  e  =  -  Bg  =  Bp  -  6^  =  (|>  -  (t.g 


11 


lOl  Q, 


When  4>^4's»  field  point  is  in  the  shadow  region;  therefore,  it 
follows  that  and  e^O  in  the  shadow  region.  The  integral  on 

the  R.H.S.  of  \12j  can  be  evaluated  via  the  method  of  steepest 
descent;  hence,  this  integral  over  the  countour  Cg  may  be  changed 
to  an  integral  along  the  steepest  descen  path  (SDP)  through  the 
saddle  point  3s  as  indicated  in  Fig.  5.  However,  in  the  process  of 
deforming  the  contour  from  Cg  to  SDP,  the  pole  at  3p  is  crossed;  it 
is  therefore  necessary  to  include  it's  residue  in  evaluating  (12). 


SDP- 

lm;3 


Re$ 


I  SHADOW  REGION 
,•  FOR  0  «f><Tr 


Fig.  5.  Location  of  Bs,  6p  and  SDP  for  the  shadow 
region  analysis. 


The  integral  in  (12)  now  becomes: 
^-ikaij*  t  -ikp[sin  B  - 

W  J 


_  ^ikp  cos  ^ 

-  -e  H(^-<}.^)  - 


^-ikoCsin  3  - 

dB  - - 

3  -  ijj 


4 


where  the  first  term  on  the  R.H.S.  of  (14)  is  the  residue  at  g  , 
and  is  the  unit  step  function:  P 

,  4,  <  (lit  region) 

,  41  >  4^  (shadow  region) 


f  0 

(15)  5  I  , 


Therefore,  the  first  term  on  the  R.H.S.  of  (14)  exactly  cancels  the 
incident  field  corresponding  to  (11)  in  the  shadow.  Approximating 

[sin  3-  in  the  exponent  of  the  integrand  on  the  R.H.S  of  (14) 
by  its  three  term  Taylor  expansion  about  the  saddle  point  gives: 


(16) 


g-ika4i 


1 


de 


g-ikp[sin  3  - 
3-4- 


-iks  -  ika  e 


2  Tii 


1 


iks[^3- 


d3 


SDP 


3  -  3. 


» 


where 

(17)  s^  i  -  a^ 


The  integral  on  the  R.H.S  of  (16)  is  now  in  a  form  which 
can  be  readily  evaluated  asymptotically  for  large  ks  by  the 
Pauli-Clemnow  modified  method  of  steepest  descent  [16,17]  for  a 
pole  close  to  a  saddle  point,  to  yield  the  following  result  for 
the  integral  in  (12). 


(18) 


m 

2 


'1 


.-ikae  „ 

+  ^  ^ - F[kL3:]  ^ 

/fnr  0  Vs 


where  the  function  F[kLa]  involves  a  Fresnel  integral  [17,18]. 


(19)  F[kLa]  =  2  i  4La  J 

/kU 


e-*’^dx  i 


in  which  one  takes  the  positive  branch  of  ^La  . 

The  distance  parameter  L  in  (19)  for  this  case  is 
(20)  L  =  s  ; 


also,  the  quantity  a  in  the  argument  kLS'  of  F  is 


It  is  noted  that  kL  is  always  the  large  parameter  in  this 
asymptotic  solution  involving  F[kLSQ.  Plots  of  the  magnitude  and 
phase  of  the  function  F[kLSG  versus  kLaf  are  illustrated  in  Fig.  6. 
The  physical  significance  of  the  geometric  quantities  6  and  s 
which  appear  in  (18)  is  evident  from  Fig.  7.  Combining  (11)  and 
(18),  one  obtains  the  following  asymptotic  result  for  I-j. 


(22a) 


Or, 


gikp  cos<t> 


gikpcos4' 


H(iJ>-4>g) 


g-ikae 

e 


F[kLa] 


>  lj)g»  0<(J)<T7. 


(22b) 


1  44" 


„-ika  6 
e _ 

6 


F[kLa] 


0<(J)<TT  . 


Turning  next  to  the  evaluation  of  in  (10b),  it  is  observed 
that  the  integrals  in  (10b)  can  be  simplified  if  the  Hy>4(kp) 

terms  appearing  therein  are  replaced  by  the  well  known  Debye 
approximation  which  is  valid  for  v(T)<kp  and  v(t:)  large.  Since 
v(t)  =  ka  +  mT  and  v(t)  0(ka)  in  the  transition  region,  the  major 
contribution  to  the  integrals  in  (10b)  occur  for  t  small.  In  this 
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d  d 

3aniiN9VM 


©  =  t  -$^>0 


Fig.  7.  Geometrical  parameters  for  the  shadow  region. 


case,  ■■X—  =  —  +  T  — —  so  that  v(t)  <  kp  is  satisfied, 

Kp  p  Kp  P 

12) 

and  the  Debye  approximation  may  be  employed  for  H  )  ukp)  as 
follows. 


H  i^!(kp)'\-  I  — 
v(Tj'  yTrkp  siny 


-ikp  sin  Y  +  iv(T)Y  +  i^ 


wi  th 


v(T)<kp,  and  cos  y  = 

K  p 


When  v(T)'^0(ka),  cos  y  ^  -  so  that  y  ^  g  ,  and 


sin  Y= 


(kp)^  -  (v(t))^ 
kp 


p 

1/2 


s' 

1/2 


kp 


—  .  Incorporating 


these  simplifications  in  the  Debye  approximation  yields: 


2  -iks  +  ij  i (ka  +  mi) 
e  ^  e 
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Next,  incorporating  (23),  and  the  relation  ()'  W^(t)  =  ^  W2(t)  +  2  i^V{T) 
into  the  integrals  in  (10b)  leads  to  the  following  expression  for  I2: 


(24a)  I2  ^  J-j^  e 


Iz  ^-iks  -  ikae  e 


.TT  - 

-  ^  / 

-  2  •' 

L  c, 


d. 


! 


^-inze 


c^+c^  ^ 


I2  can  be  further  simplified  in  a  straight  forward  manner  to  obtain: 


(24b)  I2  ^ 


e-^kae^-iks  fF  -ika  e 
“6 — F--  "’l/F  ® 


P.(C)  ^ 


where  P^(^)  is  the  Pekeris  caret  function  [9]  defined  by: 


(25)  P^(C)  H  ^ 

h 


,  ell  fijKzl 

^  _  Q  W2(t) 


and  the  parameter  5  is  defined  by 


(26)  C  =  m  e  (C>0  in  the  shadow  region). 


The  subscript  s  on  P5  is  used  to  denote  the  TM^  (or  acoustical 
soft)  case;  whereas,  the  subscript  h  on  P^  denotes  the  TEz  (or 
accoustical  hard)  case.  Values  of  P5(c)  versus  4  are  presented 

h 

in  Fig.  8  in  terms  of  related  functions  p*(c)  and  q*(c). 
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e 


Fig.  8(b).  Plot  of  versus  5  based  on  Logan's 

tabulated  data  [9]  for  q(5). 
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Combining  the  results  of  (22b)  and  (24b),  one  obtains  the 
following  asymptotic  expression  for  U  (of  Eq.  (9))  which  is 
valid  at  the  observation  point  in  the  shadowed  portion  of  the 
transition  region. 


(27) 


U(P^)'v. 


/2nk 


-ikae 


0 


-  F[kLa] 


-  m  Ps(c)  ;  e  >  0,  4>  > C>0. 

h 


The  above  result  for  U  may  be  interpreted  geometrically  as 
a  field  which  after  being  launched  by  the  incident  wave  at  Q] 
propagates  along  the  path  QiQ2Ps  to  the  field  point  Pg,  as  indicated 
in  Fig.  7.  Q]  and  are  tangent  points  which  correspond  to  the 
point  of  grazing  incidence  on  the  cylinder,  and  the  point  of  tan¬ 
gential  shedding  of  the  diffracted  wave  from  the  surface,  respec¬ 
tively. 


One  may  now  examine  the  result  of  directly  extending  the  analysis 

presented  above  to  lit  region.  In  the  lit  region,  and  as 

s  P  s 

indicated  in  Fig.  9  (note:  Fig.  5  shows  in  shadow  region). 

Consequently,  H(<f)-<^2)  =  0;  also,  e<0  (^<0)  fori}><$^,  and  the  field  U 
at  the  point  Pj^  in  the  lit  region  becomes  via  (22a)  and  (24b): 


(28) 


U(P^)  -v  e^*^^ 


-  F[kLa] 


j-ikao 


h 


e 


-iks 


/s 


e<0;  ?<0. 
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Fig.  9.  Location  of  Bp  and  SDP  for  the 
lit  region. 


Equations  (27)  and  (28)  may  be  represented  compactly  in  a 
single  expression  as  follows: 


(29) 


U(P)  -v  H(-0) 


/TTk 


-  m  ||  P  (C) 

h 

2 

with:  kLa  =  kL'^  ,  €  =  me  ,  L  =  s  ;  also, 

P  =  P|^  ,  for  5<0,  e<0,  and  'l><ii>g(in  the  lit  region), 

P  =  Pg  ,  for  ^>0,  0>O,  and  (in  the  shadow  region), 
r 

1 ,  0<C  for  P  =  P|^ 

0,  e>0  for  P  =  P„ 


H(-e)  = 


The  first  term  on  the  RHS  of  (28)  and  (29)  is  just  the  incident 
field  uHP).  The  geometrical  interpretation  for  the  propagation 
of  the  field  L)  in  the  lit  region  is  presented  in  Fig.  10;  this 
interpretation  is  simply  an  extension  of  that  in  Fig.  7  for  the 
shadow  region;  more  will  be  said  about  this  later. 


0  -  I//  -  ^3  <  0 


Fig.  10.  Geometrical  parameters  for  the  lit  region. 


The  result  in  (29)  will  be  shown  to  be  finite  and  continuous 
across  SBi .  As  one  approaches  SBi,  and  e-+0;  hence, 

2 

a  =  ^"^0  3nd  cX).  The  limiting  forms  of  F[kLal  and  as  e, 

are  [18,19]:  ^ 


irkLa 


-  2 


kLa 


'r 


.  /  7T 

^(4 


kLa) 


f;-Kl 


(30a)  F[kL^] 


e 


for  a  given  kL  as  a^O;  and 


2  /rr^ 

+  (small  contribution  from  p*  or  q*  which  is 
continuous  w.r.t.  ^),  as  ^-*0. 

■H  -TI 

“^T  “It 

(Note:  the  p*(i;)e  and  q*(c)e  functions  are  plotted  in  Fig.  8.) 

Employing  the  above  limiting  forms  of  (30a;  30b)  in  (29)  as 
e,  ^'*■0  from  both,  the  lit  (e,  ?<b)  and  shadow  (e,  c>0)  directions, 
it  is  easily  verified  that  U  in  (29)  is  finite  and  continuous 
across  SB]*,  in  fact,  exactly  at  SB-j ,  L)(P)  becomes: 


(30b)  P^(0  = 


JP*(5)1 
q*(?) I 


2  /ttC 


(31) 


U(P) 


m 


SB, 


for '! 


TM^  case  ! 


j^TE^  case 


1’ 


'P*(0)1  pikx 
q*(0)  j  j/xl 


where  x  =  p  cos  ijig  =  -  s  along  S8i  (a  similar  result  holds  for  SB2). 

The  first  term  on  the  R.H.S.  of  (31)  is  simply  one-half  the  incident 
field  on  the  shadow  boundary  as  in  the  diffraction  by  a  half-plane; 
the  second  term  is  dependent  on  surface  properties. 

Turning  now  to  the  evaluation  of  the  field  in  the  shadow  region, 
but  away  from  SBi  and  exterior  to  the  transition  region  (i.e.,  in 
region  III  of  Fig.  1),  one  notes  that  FEkLa]"^!  since  e,  and  hence 
kLat  become  large  in  this  case  (see  Fig.  6);  only  the  term  containing 
p5(^)  in  (29)  now  remains  significant.  Far  from  SB-j  the  integral 

h  « 

for  Pg(c)  in  (25)  may  be  replaced  by  a  rapidly  convergent  residue 
h 

series  [9]: 


23 


(32) 


P3(C) 


.  TT 

I  e 


i.TT 


„  -  - — TT  ,for  $=1  or  TM,  case 

^  n  2[Ai'(-qJ]2  " 


=  < 


.  5tt 


5>0 


-'T 


2  £  ® 


— z - r=r  *  r“  or  TE 

n  2[Ai(-q)]2q„  " 


where  Ai  is  the  Mi  Her- type  Airy  function*  (which  is  related  to  the 
Fock-type  Airy  function  [9]),  and  Ai'_is  it's  derivative  with 
respect  t£  the  argument.  The  q^  and  qn  are  roots  of  Ai(-qp)  =  0 
and  Ai'(-qn)  =  0,  respectively  for  n  =  1,2,3,  These  roots 

are  tabulated  in  [9]. 


When  the  above  residue  series  representation  of  (32)  is 
employed  in  P^(c)  of  (29)  for  kL3f>10  (i.e.,  as  F[kLa>l),  it 

h 

can  be  readily  shown  that  (29)  reduces  to 


(33) 


U(P,) 


u'(Q,) 


)  e 


-ikae  - 


a^as  ? 


^-iks 

TT 


s  s 

where  and  are  the  Keller  diffraction  and  attenuation  coef¬ 
ficients  [3]  associated  with  the  GTD  surface  ray  modes,  respec- 
s  s 

tively,  with  d||(Qi)  =  ^^(^2^  circular  cylinder.  The 

s  s 

h  h 

superscripts  in  D"  and  a"  correspond  exactly  to  the  subscripts  in 
^  n  n 

P3(^)  for  the  two  polarizations;  also,  the  subscript  n  denotes  the 
h 


“  2 

Ai(T)  =  y{T)//n  =  g~i(Tt+t  /3)  .  = 


dAi (t  ) 

3t 


case 
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surface  ray  modal  index.  The  U^(Q, )  in  (33)  represents  the  incident 
field  at  Q]  and  is  unity  for  this  problem.  Thus,  in  the  deep  shadow 
(i.e.,  corresDondinq  to  region  III  of  Fiq.  1).  (29)  reduces  to  thp 
GTD  surface  diffracted  ray  field  given  by  (33)  as  it  should.  The 
geometrical  interpretation  of  the  GTD  field  in  (33)  is,  of  course, 
identical  to  that  in  Fig.  7.  For  the  sake  of  convenience,  and 
completeness,  the  diffraction  and  attenuation  coefficients  in  (33) 
are  listed  below. 


,2  /7  2-5/^  1/3^-^./12 


Ai'(-qn) 


.2  2-5/^  1/3  ^-W12 


k’'®  q, 


h  .  ^ 

n  p 


%  ^iTT/ef^p 


B.  Transition  Region  -  Iluminated  Zone 

Although  the  result  in  (29)  is  applicable  to  both  the  shadow 
and  lit  portions  of  the  transition  region,  it  does  not  reduce  to 
the  geometrical  optics  field  in  the  lit  region.  Consequently,  the 
result  in  (29)  is  expected  to  be  accurate  in  the  shadow  region 
both  near  and  far  from  SB] ;  whereas,  it  is  expected  to  be  accurate 
only  for  e  small  in  the  lit  region  U,  e<0).  This  deficiency  of 
(29)  for  the  lit  region  where  eis  not  small  is  briefly  discussed, 
and  an  improved  solution  which  overcomes  this  difficulty  is  presented. 

For  the  shadowed  part  of  the  transition  region,  the  result  in 
(29)  has  the  geometrical  interpretation  illustrated  in  Fig.  7 
which  is  identical  to  the  GTD  surface  diffracted  ray  picture  as 
shown  in  Fig.  1;  consequently,  it  is  not  too  surprising  that 
Eq.  (29)  does  indeed  reduce  exactly  to  the  GTD  solution  of  Eq. 

(33)  in  the  deep  shadow.  On  the  other  hand,  the  geometrical 
interpretation  of  the  scattered  field  in  (29)  for  the  lit  region 
(0<O)  is  illustrated  in  Fig.  10  which  constitutes  a  direct 
generalization  of  the  interpretation  in  Fin.  7  (for  9>0);  clearly, 
the  ray  geometry  in  Fig.  10  is  in  disagreement  with  the  geometrical 
optics  reflected  ray  picture  of  Fig.  11  which  is  known  to  be  valid 
in  the  lit  region.  Consequently,  one  might  anticipate  the  failure 
of  Eq.  (29)  to  reduce  to  the  geometrical  optics  field  in  the  lit 
region  as  a  result  of  the  differences  between  the  scattered  field 
propagation  mechanisms  of  Figs.  10  and  11. 
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Fig.  11.  Reflected  ray  system. 


From  Fig.  10,  it  is  observed  that  the  scattered  field  at  P| 
traverses  the  path  Q1Q3PL  after  being  launched  by  the  incident 
wave  at  Q-i.  This  interpretation  requires  the  field  to  propagate 
backward  from  Q-[  to  Qo  and  then  shed  to  Pl  along  the  backward 
tangent  at  Q3;  in  contrast,  the  geometric  optical  reflected  ray 
trajectory  touches  the  surface  only  at  Qr  as  in  Fig.  11.  The 
paths  in  Figs.  7  and  11  satisfy  the  extended  Fermat's  principle; 
whereas,  the  ray  path  in  Fig.  10  does  not.  The  inability  of  Eq. 
(29)  to  reduce  to  the  geometrical  optics  field  may  be  explained 
on  the  basis  of  the  approximation  in  (23)  for  the  Hankel  function 

Hy|^|(ko)  which  is  employed  in  the  evaluation  of  the  transition 

region  field;  more  will  be  said  about  this  shortly. 

For  the  sake  of  completeness,  the  behavior  of  Eq.  (29)  in 
the  deep  lit  region  is  examined.  As  before,  F[kL&]-^l  as  the 
field  point  moves  away  from  SB-j  in  bqth  the  lit  and  shadow  zones, 
thereby  leaving  the  term  containing  Pg(c)  as  the  dominant  con- 

h 

tributor  to  the  scattered  field.  Hence,  in  the  lit  region  far 
from  SB]  where  ^«0,  Eq.  (29)  reduces  to 


(34) 


U(Pl)  ^  H(-e) 


-iks 


"Jr 

h  ^ 


e<0,  c«0. 


When  C«0»  PgCc)  has  the  following  asymptotic  value  [9] 


(35)  ^5(0  [1  +  0(r^)]  ;  C«0. 

h 

Incorporating  (35)  into  (34),  one  arrives  at  the  following 


(36)  U(Pl)  u'(Pl)  + 


ikaCje]  - 


/s 


'TM. 


fori  ^  vcase  when  5«0. 


On  the  other  hand,  the  geometrical  optics  field  represents 
the  correct  asymptotic  result  far  from  SB]  in  the  lit  region, 
and  its  field  associated  with  the  incident  and  reflected  rays 
is  well  known  to  be 


(37)  U(Pl)  'V  u’(Pl)  + 


-ikA 


fPr 


ika  cos  e 


where  U^Pl)  =  e^'^^  *  and  U^Qr)  =  e""“  (see  Fig.  11). 

The  first  term  on  the  RHS  of  Eq.  (37)  is  the  incident  field  at  ?[_, 
and  the  second  term  denotes  the  reflected  ray  field.  Qr  is  the 
point  of  reflection  on. the  surface  as  shown  in  Fig.  11,  and  the 
{+}  sign  in  front  of  U''(Qr)  denotes  the  value  of  the  surface 


reflection  coefficient  at  Qr  for  the 


distance  t 
associated 


IS  given 
with  the 


by  4  =  q 
reflecte 


5%i 


whereas , 
ield  is 
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case.  The  spatial 


the  caustic  distance  0, 


(38) 


a  cos  0^ 


The  quantities  i  and  are  also  shown  in  Fiq.  11. 

One  notes  that  Eq.  (36)  and  Eq.  (37)  are  nearly  equal  only  for 
0  very  small.  In  general,  therefore,  the  result  of  Eq.  (29)  fails 
to  reduce  the  geometrical  optics  result  when  c<<0  (this  failure  is 
similar  to  that  which  occurs  in  the  results  of  [7],  [8],  and  [14]); 
consequently,  in  order  to  overcome  this  limitation,  one  must  employ 

(2) 

an  approximation  for  lit  region  which  is  different 

from  that  in  (23),  and  which  leads  to  a  geometrical  interpretation 
for  the  scattered  field  in  the  lit  region  which  is  in  agreement 
with  the  reflected  ray  geometry  of  Fig.  11. 

The  approximation  of  Hyf|](kp)  in  (23)  is  valid  for  v(t) 0(ka) ; 
the  condition  v(T)'^0(ka)  holds  in  the  transition  region  near  SBi 
and  also  in  the  deep  shadow  zone  where  only  the  first  couple 
of  GTD  surface  ray  modes  (corresponding  to  n  =  1,  2,  etc.,  in  (32) 
and  (33))  are  necessary  to  describe  the  field  accurately;  whereas, 
deep  in  the  lit  region  v(t)  can  be  significantly  different  from  ka. 
Thus,  Eq.  (23)  leads  to  the  result  of  Eq.  (29)  which  is  accurate  in 
the  shadow  region,  but  it  is  not  expected  to  be  accurate  in  the  lit 
region  if  the  field  point  is  not  close  to  SBi. 

The  Fock-type  Airy  terms  Q  V(t),  and  Q  W^(t)  in  (10b)  which 

(2)  ^ 

approximate  Jy(ka)  and  Hy'  '(ka)  for  v'^'O(ka)  reduce  uniformly  to 
the  Debye  asymptotic  approximations  when  x  becomes  large  (x  = 

—  )  ;  hence,  in  contrast  to  the  approximation  for  Hy(^)('^p) 

in  (23),  these  Fock-type  Airy  functions  remain  valid  as  v(x)  becomes 
different  from  ka  (i.e.,  for  jv  -  kal>  0(m))  deep  in  the  lit  region. 
Instead  of  the  approximation  in  (23),  one  could  also  approximate 
(2) 

^v(x)(kp)  Fock-type  Airy  function 

W2  as  follows,  provided  k[p  -  a]«ka. 

(39)  ^  '^2  '»  l<[p-a]^<l<a. 

^  m  /tt 


But,  employing  the  near  zone  approximation  of  Eq.  (39)  in  I2  of 
Eq.  (10b)  leads  to  integrals  which  are  complicated  functions  of 
two  parameters  (namely,  and  [p  -  a]);  on  the  other  hand,  the 
result  in  (29)  is  in  terms  of  simpler  F  and  integrals  which  are 

h 

each  functions  of  only  a  single  parameter. 
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A  representation  for  U  in  the  lit  region  which  is  interpretable  in 
terms  of  the  desired  reflected  ray  geometry  of  Fig.  11  and  which  also 
is  in  terms  of  the  F  and  functions  may  be  obtained  in  the  follow- 

^  f  21 

ing  manner  by  incorporating  an  approximation  for  which  unlike 

(23)  is  valid  even  if  v  becomes  quite  different  from  ka  in  the  lit  region. 
One  begins  by  introducing  the  Debye  approximation  for  Hy|^|(kp)  which 
was  previously  indicated  as 


(40) 


i  (2) 

v(t) 


(kp) 


Tfkp  sin  Y 


-ikp  sin  Y+  iv(T)Y  +  i 
e 


IT 

T 


v(t)  =  kp  cos  Y  and  v(t)<  kp. 


The  Condition  that  v(T)<kp  be  satisfied  in  the  lit  zone  is  verified 
in  the  Appendix  where  the  geometric  optical  reflected  field  result 
(corresponding  to  the  second  term  on  the  R.H.S.  of  Eq.  (37))  is 
developed  in  detail.  Referring  to  this  development  in  the  Appendix, 
it  is  seen  that  in  the  lit  region  where  v<ka,  the  following  approxi¬ 
mations  are  valid: 


(41a)  v(t)  ^  Vq  =  ka  cos  =  kp  cos  y^  ;  Y  ^  Yp 

(Alb)  sin  Yq  =  cos  e\  since 

(41c)  p  sin  Y  ^  p  sin  Yq  =  *'■•■  a  sin  y^  =  c  +  a  cos 

(41d)  ^0  "  2  Yq 

(Note:  the  geometrical  signficances  of  (41)  is  shown  in  Fig.  A-I 
in  the  Appendix. ) 

The  quantity  ip  has  already  be  defined,  and  y^  |s  introduced 
in  the  Appendix.  The  subscript  o  in  Vo(t),  yq  and  yq  refers  to 
the  geometrical  optics  reflected  field  saddle  point  condition 
which  is  also  discussed  in  the, Appendix.  Incorporating  v(t)y  ^ 
v(t)  yqJ  p  sin  Y  ^  +  a  cos  ei  for  the  phase;  and  p  sin  y%i  for 
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in  the  exponent.  This  approximation  is  valid  in  the  transition 
region.  One  notes  that  the  approximation  in  (45a)  and  (45b)  follow 
directly  from  the  expansion 


(45c) 


■Y 


=  Sin 


•1  =  J_ 

2m  2m 


6(2m)‘ 


40 (2m)' 


wr 


with 


<1. 


Employing  (45a)  and  (45b)  into  (43)  yields 


(46) 


T  .  (2  Q-Ikt  -ika  cos  „  „i2ka  cos  eV 

lo  -  -  e  me 


■i(5')^/12 


-I: 


e -  +  P  (5') 

2^C-  ^ 


which  is  valid  in  the  far  zone.  In  the  lit  region  far  from  SB], 
C'<<0,  and  ^^(C)  may  be  replaced  by  (35)  to  yield 

h  •  I — 

(47)  ij  .  5  ®  ,1^  e-*"'  , 


which  does  indeed  agree  with  the  far  zone  geometrical  optics 

reflected  field.  The  term  e  corresponds  exactly  to 

d^(0^)>  and  {+}denotes  the  reflection  coefficient  values  as  in  (37). 

It  is  now  evident  that  one  can  employ  ray  optics  to  generalize  (46) 
so  that  it  yields  the  near  zone  geometrical  optics  reflected  field 
as  in  (37).  According  to  ray  optics,  the  far  zone  energy  spread 

Fr 

factor  Jj—  in  (47)  must  be  simply  replaced  by  its  corresponding  near 


zone  factor 


pz 

VPf,  + 


(as  in  (37))  for  the  reflected  ray  tube  in  order 


to  extend  (47)  to  the  near  zone  case;  the  latter  is  tantamount  to 

replacing  —  in  (46)  with  —  -  -  thereby  generalizing  (46)  to  the 

/I  /p“+T 

near  zone  case. 
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To  be  consistent,  it  is  also  necessary  to  modify  the  result 
for  Ii  of  (10a)  for  the  lit  renion  in  the  same  manner  as  done  above 
for  I2  of  (10b).  Without  giving  all  the  details  (since  these  involve 
essentially  the  same  physical  arguments  as  in  the  treatment  of  I2), 
the  result  for  I]  in  the  lit  region  becomes 


with 


(48b)  L'  t  ,  for  C'<0 


Finally,  combining  (46)  (together  with  it's  near  zone  modifica¬ 
tion)  and  (48),  one  obtains  the  required  result  for  U  in  the  lit 
region;  namely. 


(49)  U(Pl)  ^  U^P^)  +  U^O^) 


-n  e  g 

/zTk  c' 


m3 


•  J  1  -  F[kL' 


2m  J 


-ikt 


/i.+ 


-  m 


/! 


-i(C')' 


P,U') 


„-ik£. 
e _ 

/t+  n 


;  s'<o. 


It  is  easily  verified  that  far  from  SB^  (5'<<0),  F[kL' 


]-l 


and  P^(^')  behaves  as  in  (35)  so  that  in  the  deep  lit  regions  Eg. 
h 

(49)  reduces  exactly  to  the  geometrical  optics  result  of  Eq.  (37). 
Furthermore,  Eg.  (49)  is  also  valid  for  9'-^/2  (i.e.,  near  SBi) 
since  it  agrees  in  this  limiting  case  with  the  expression  in  (29) 
for  ^.<0.  In  addition,  Eq.  (49)  reduces  to  Eq.  (31),  thereby 
assuring  the  boundedness  and  continuity  of  the  field  across  SB], 
Thus,  the  result  for  U  in  (49)  for  the  lit  region  reduces  to  the 
geometrical  optics  field  deep  in  the  illuminated  region  where  the 
latter  is  indeed  valid,  and  at  the  same  tine  it  has  the  proper 
behavior  required  to  ensure  the  continuity  of  the  field  across 


III.  NUMERICAL  RESULTS 

In  this  section,  numerical  results  are  presented  for  the  near 
zone  total  field  surrounding  a  perfectly-conducting  circular 
cylinder  excited  by  a  plane  wave.  The  GTD  method  is  employed  to 
describe  the  field  exterior  to  the  cylinder  in  terms  of  the 
geometrical  optics  incident  and  reflected  rays  for  the  lit  region, 
and  in  terms  of  the  surface  diffracted  rays  for  the  shadow  region, 
respectively;  these  results  for  the  TE2  (hard  boundary)  case  are 
indicated  in  Figs.  12(a)  -  20(a)  for  different  size  cylinders 
and  near  field  distances.  It  is  seen  that  the  GTD  solution  shown 
by  a  dashed  line  becomes  discontinuous  at  the  shadow  boundaries. 

On  the  other  hand,  the  present  transition  region  solution  (based 
on  (27)  for  the  shadow,  and  on  (49)  for  the  lit  region)  as  indicated  by 
a  solid  line  is  continuous  across  the  shadow  boundaries,  and  it 
blends  smoothly  into  the  GTD  solution  outside  the  transition  region. 
These  GTD  and  transition  region  solutions  corresponding  to  Figs. 

12(a)  -  20(a)  are  combined  into  a  single,  composite  GTD-transition 
solution  which  is  shown  by  a  smooth  solid  curve  in  Figs.  12(b)  - 
20(b),  respectively.  The  composite  GTD-transition  results  are 
compared  against  results  which  are  based  on  an  exact  eigenfunction 
series  solution  for  these  cases;  the  eigenfunction  results  are 
shown  by  dots  in  Figs.  12(b)  -  20(b)  for  these  cases.  Figures 
12(b)  -  20(b)  indicate  that  the  comparison  between  the  numerical 
results  obtained  from  the  exact  eigenfunction  solution,  and  the 
composite  GTD-transition  solution  is  very  good,  thereby  confirming 
the  accuracy  of  the  uniform  result  for  the  transition  region 
presented  in  section  II.  Additional  comparisons  of  the  composite 
GTD-transition  solutions  with  the  corresponding  eigenfunction 
solutions  are  indicated  in  Figs.  21-22  for  the  TM^  (soft  bound¬ 
ary)  case,  respectively;  again  the  agreement  is  very  good  for  these 
cases.  It  is  seen  that  the  accuracy  of  the  transition  region  solu¬ 
tion  for  very  small  ka  and  kp  is  as  good  as  the  accuracy  of  the 
geometrical  optics  result.*  As  ka  and  kp  increase,  our  asymptotic 
transition  field  solution  becomes  increasingly  accurate  as  do  the 
GTD  (surface  diffracted  and  geometric  optical)  ray  fields  which 
also  are  asymptotic  approximations.  For  calculating  the  GTD  field 
in  the  shadow  region,  the  inclusion  of  only  two  surface  ray  modes 
are  found  to  be  sufficiently  accurate  (see  (33)).  Unfortunately, 
there  appears  to  be  a  small  error  in  the  p*  function  plotted  in 
Fig.  8(a);  this  error  is  noticeable  for  ranges  of  C  which  occur  for 
small  kp  and  ka  in  the  TM7  case.  Presently,  improved  curves  for  the 
amplitude  and  phase  of  p*(^)  are  being  obtained. 


*In  the  GTD  calculation  for  the  lit  region  (or  the  geometrical  optics 
region),  the  surface  rays  that  have  crept  around  the  opposite  side 
of  the  cylinder  have  been  neglected;  but  the  field  of  these  rays  is 
included  with  the  transition  field  in  the  calculation  shown  by  the 
solid  line. 
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TRANSITION 

FUNCTION 


p  =0.5X 
HARD  b.c. 


Fig.  12(a),  Field  surrounding  a  circular  cylinder  illumi¬ 
nated  by  a  TE^  plane  wave.  Calculations  based 
on  GTD  and  the  transition  solution. 
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a  =0.4X 
p  =0.5X 
HARD  b.c. 


Fig,  12(b).  Comparison  of  the  composite  GTD  plus  transition 
solution  of  Fig.  12(a)  with  the  exact  solution. 
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13(a).  Field  surrounding  a  circular  cylinder  illumi¬ 
nated  by  a  TE2  plane  wave.  Calculations  based 
on  GTD  and  the  transition  solution. 


HARD  b.c. 


Fig.  14(a).  Field  surrounding  a  circular  cylinder  illumi¬ 
nated  by  a  TE^  plane  wave.  Calculations  based 
on  GTD  and  the  transition  solution. 
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HARD  b.c. 


Fig.  15(a).  Field  surrounding  a  circular  cylinder  illumi¬ 
nated  by  a  TE^  plane  wave.  Calculations  based 
on  GTD  and  the  transition  solution. 
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Fig.  16(a).  Field  surrounding  a  circular  cylinder  illumi¬ 
nated  by  a  TE^  plane  wave.  Calculations  based 
on  GTD  and  the  transition  solution. 
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270^  yo  =  2  X 

HARD  b.c. 

Fig.  17(a).  Field  surrounding  a  circular  cylinder  illumi¬ 
nated  by  a  TE^  plane  wave.  Calculations  based 
on  GTD  and  the  transition  solution. 
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GTD  +TRANSITION 
SOLUTION 


HARD  b.c. 


Fig.  17(b).  Comparison  of  ,  o  composite  GTD  plus  transition 
solution  of  Fig.  17(a)  with  the  exact  solution. 
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HARD  b.c. 

18(a).  Field  surroundinp  a  circular  cylinder  illumi¬ 
nated  by  a  TE2  plane  wave.  Calculations  based 
on  GTD  and  the  transition  solution. 
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GTD  +TRANS1TI0N 
SOLUTION 


p  =3X 
HARD  b.c. 


Fig.  18(b).  Comparison  of  the  composite  GTD  plus  transition 
solution  of  Fig.  18(a)  with  the  exact  solution. 
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Fiq.  19(d).  Field  surrounding  a  circular  cylinder  illumi¬ 
nated  by  a  TE^  plane  wave.  Calculations  based 
on  GTD  and  the  transition  solution. 
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Fig,  20(a).  Field  surrounding  a  circular  cylinder  illumi¬ 
nated  by  a  TEz  plane  wave.  Calculations  based 
on  GTD  and  the  transition  solution. 
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a  =  I.59X 
p  =  4.76  X 
SOFT  b.c. 


Fig.  21(a).  Field  surrounding  a  circular  cylinder  illumi¬ 
nated  by  a  plane  wave.  Calculations  based 
on  GTD  and  the  transition  solution. 
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p  =7.95X 
SOFT  b.c. 


Fig.  22(a).  Field  surrounding  a  circular  cylinder  illumi¬ 
nated  by  a  TM^  plane  wave.  Calculations  based 
on  GTD  and  the  transition  solution. 


SOFT  b.c. 


Fig.  22(b),  Comparison  of  the  composite  GTD  plus  transition 
solution  of  Fig.  22(a)  with  the  exact  solution. 


IV.  SUMMARY  AND  DISCUSSION 

A  uniform  asymptotic  result  is  developed  for  describing  the 
fields  in  the  transition  regions  adjacent  to  the  shadow  boundaries 
which  are  associated  with  the  problem  of  plane  wave  diffraction  by 
a  perfectly-conducting,  circular  cylinder.  This  result  yields 
a  finite,  continuous  solution  for  the  total  field  across  the  shadow 
boundaries,  and  exterior  to  the  transition  regions  it  reduces  to 
the  usual  GTD  solution  (in  terms  of  the  surface  diffracted  fields 
and  geometrical  optical  ray  fields)  for  the  shadow  and  lit  zones, 
respectively. 

Sunmarizing  this  result  for  convenience,  one  expresses  the 
total  field  U  exterior  to  the  cylinder  via  (27)  and  (49)  as 

TRANSITION  REGION  -  SHADOW  ZONE 

—  •  ‘H’ 


with  c  =  me>0. 


(51a)  U(P^)  U'(Q^) 


2/i75  'L 


1  -  F 


where  t  =  as  = 


ka  i 


1/3 


The 


Q^Q2,  s  =  Q2P5  as  in  Fig.  7;  and  m  = 

functions  F[ks  and  P  (^)  are  tabulated  in  Figs.  6  and  8, 

2-^  h 

respectively;  one  should  note  that  Ps(c}  is  associated  with  the 
acoustic  soft  boundary  case  (TM^),  and  Ph(c)  is. associated  with 
the  acoustic  hard  boundary  case  (TEz).  Also,  UMO-j)  denotes  the 
value  of  the  incident  field  at  Q-]. 


TRANSITION  REGION  -  LIT  ZONE 


(51b)  U(P|_)  ^  U^P^)  +  U^Qr) 
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2  ^  ll^ 


\h 

with  5'  =  -2m  cos  6^  <0  , 


kn 


,2' 


2m 


|pr  ^ 


-ikt 


i  2  ^ 

where  ^  =  OrP^,  ^  ;  and  m,  F[kA  ,  Pg(5')  , 

2m 


56 


and  u’  have  the  same  meaning  as  before.  The  geometrical  quantities 
)t,ei  and  the  point  of  reflection  Qr  are  indicated  in  Fig.  11.  One 

^  i— 

notes  that  +  1_  =  |P* jj,  e  in  (51a)  and  (51b), 

respectively.  The  functions  p*  and  q*  are  plotted  in  Fig.  8. 

Equation  (51)  is  valid  for  both,  the  near  and  the  far  zone; 

however,  in  the  near  zone  case,  the  field  point  must  lie  exterior 
to  the  surface  boundary  layer.  Within  the  surface  boundary  layer, 

(51a)  tends  to  become  singular  (as  s-'-O)  in  exactly  the  same  manner 

as  the  GTD  surface  diffracted  ray  field  of  (33)  which  becomes 
singular  as  s-^0.  In  order  to  remove  this  limitation,  one  must 
employ  (39)  instead  of  (23)  in  the  analysis  as  s->-0;  however,  the 
case  s->0  is  not  considered  in  this  paper.  One  could  employ  the 
results  of  Ivanov  [11]  which  allows  the  field  to  be  calculated  near 
the  surface  (within  the  surface  boundary  layer)  as  s-^0.  However, 
as  mentioned  earlier,  his  result  [11]  is  in  terms  of  the  canonical 
Fock  integral  which  is  complicated  for  numerical  calculations  except 
at  s=0;  when  s=0,  the  canonical  Fock  integral  reduces  to  the  single 
parameter  Fock  current  integral  which  is  well  tabulated.  It  would 
be  worth  investigating  an  alternative  approximate  solution  which 
matches  our  solution  of  (51)  outside  the  surface  boundary  layer, 
but  which  reduces  to  the  Fock  currents  on  the  surface  (s=0),  and 
at  the  same  time  is  such  that  it  is  simple  for  numerical  processing. 

The  result  in  (51)  yields  a  finite  and  continuous  field  solution 
across  the  shadow  boundaries.  Away  from  the  shadow  boundaries,  and 
exterior  to  the  transition  region  (51a)  reduces  to  the  GTD  surface 
diffracted  ray  field  of  (33);  whereas,  (51b)  reduces  to  the  geometrical 
optics  incident  and  reflected  ray  fields  of  (37).  The  F[*]  term  in 
(51)  dominates  in  the  region  at  and  near  the^ shadow  boundaries; 
whereas,  far  from  the  shadow  boundaries  the  P^{’)  term  dominates. 

Since  the  F(*)  and  P^(*)  functions  are  well  tabulated  [9],  the  result 

in  (51)  is  simple  to^use  which  makes  it  well  suited  to  solve  practical 
problems.  Furthermore,  this  result  can  be  readily  generalized  to 
convex  surfaces  of  variable  curvature  in  both  two  and  three  dimensions 
by  employing  the  local  nature  of  high  frequency  fields  and  recipro¬ 
city.  These  extensions  will  be  presented  in  a  separate  paper. 

Finally,  the  field  at  (p,<{>)  due  to  TM^  (or  TE^)  plane  wave  scattering 
by  the  cylinder  is  directly  related  via  reciprocity  to  the  problem 
of  calculating  the  radiation  pattern  (along  the  incident  field  direc¬ 
tion)  of  an  electric  (or  magnetic)  line  source  at  (p,(j))  in  the  presence 
of  the  same  cylinder;  the  latter  problem  is  of  interest  for  estimating 
the  radiation  patterns  of  antennas  mounted  off  the  fuselage  of  a 
typical  aircraft. 
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APPENDIX 


The  result  for  the  geometrical  optics  field  reflected  from 
a  perfectly-conducting  cylinder  which  is  illuminated  with  a  plane 
wave  is  developed  in  this  Appendix;  this  development  is  useful  to 
the  analysis  of  section  II  B.  Although  this  development  is  avail¬ 
able  in  many  texts  (e.g.,  [15^,  it  is  rederived  here  in  terms  of 
the  notation  employed  in  section  II  so  that  the  connection  between 
the  results  derived  here  with  the  analysis  in  section  II  B  may  be 
established  in  a  convenient  fashion.  Only  the  TM^  case  is  treated 
here;  the  results  for  the  TE^  case  are  exactly  the  same  except  for 
a  change  of  sign. 

Referring  to  Fig.  2,  it  is  seen  that  the  total  field  exterior 
to  the  cylinder  in  this  case  is  (see  (6b)): 
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For  the  TM^  case,  Q  =  1 ,  and  the  first  term  on  the  RHS  of  (Al) 
yields  the  incident  field  [15];  on  the  other  hand,  the  second  term 
on  the  RHS  of  (Al )  yields  the  reflected  field  for  v<ka.  One  is 
concerned  only  with  this  second  term  which  is  denoted  below  by 


(A2) 
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Employing  the  Debye  approximations  in  the  Hankel  functions 
leads  to: 


(A3) 


where 

(A4) 
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V  =  ka  cos  Y  =  kp  cos  y  » 
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and  the  contour  of  integration,  C^,  is  indicated  below  in  Fig  A-I 
Also,  rp  =  (I4.I  -  Tr/2)  as  before. 
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The  following  relations  hold  as  a  consequence: 


(A8) 

(A9) 


The 

in  Fig.  A 


1  .  . 
cos  e  =  sin  Yq  ,  when  v  = 


p  sin  y  /  *  =  p  sin  Yq  =  ^  +  a  sin  Yq  =  +  a  cos  e’ 


Fig.  A-II.  Reflected  ray  geometry. 


geometrical  significance  of  (A7),  (A8),  and  (A9)  is  shown 
■II.  Also, 

f(v^)  =  a  cos  e’  -  ^ 

f'(v„)  .0 

f..(  )  .  2  a  cos  e’' -  _ 

(ka  cos  e  )(kp  sin  Yq) 


where  the  primes  denote  differentiation  w.r.t.  v 


Incorporating  the  above  expressions  into  the  saddle  point 
result  for  (A5)  gives 
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One  finally  obtains 
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- -  e  ,  for  TM^  case 


which  in  the  notation  of  section  II  B  reduces  to 


(A12)  .  u'(q„)  ;  R  =•?!; 

u  V  r  . 


which  is  the  required  reflected  geometrical  optics  field. 


